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Abstract 

Two main seismic features characterize the Earth's inner core: a North-South polar anisotropy and an East- West 
asymmetry of P-wave velocity and attenuation. Anisotropy is expected if shear deformation is induced by convective 
motions. Translation has recently been put forward as an important mode of convection of the inner core. Combined 
with a simple diffusive grain growth model, this mechanism is able to explain the observed seismic asymmetry, but not 
the bulk anisotropy. The source of anisotropy has therefore to be sought in the shear motions caused by higher modes 
of convection. Using a hybrid finite-difference spherical harmonics Navier-Stokes solver, we investigate the interplay 
between translation and convection in a 3D spherical model with permeable boundary conditions at the inner core 
boundary. Three parameters act independently: viscosity, internal heating and convection velocity in the outer core. 
Our numerical simulations show the dominance of pure translation for viscosities of the inner core higher than 10 20 Pas. 
Translation is almost completely hampered by convective motions for viscosities lower than 10 18 Pas. Between these 
values, translation and convection develop, but convective downwellings are restricted to the coldest hemisphere where 
crystallization occurs. On the opposite side, shear is almost absent, thereby allowing grain growth. We propose that the 
coexistence of translation and convection observed in our numerical models leads to a seismic asymmetry but localizes 
deformation only in one hemisphere. 
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Introduction 

The image of the inner core growing slowly at the cen- 
ter of the Earth by gradual cooling and solidification of the 
surrounding liquid outer core is being replaced by the more 
vigorous image of a "deep foundry" (Buffctt, 2011), where 
melting and crystallization rates exceed by many times the 
net growth rate (Monnereau et al., 2010; Alboussiere et al., 
2010; Gubbins et al., 2011) 

During seventy years, the analysis of compressive waves 
(P waves) and free oscillations excited after large earth- 
quakes have been depicting a more and more complex 
structure of the inner core. It appears anisotropic, with 
a fast axis parallel to Earth's spin axis. It is also asym- 
metric: within the outermost 100 km, velocity and at- 
tenuation of P waves increase from the hemisphere facing 
America (West) to the one facing Asia (East) (Tanaka and 
Hamaguchi, 1997) ; anisotropy also seems stronger in the 
Western hemisphere than in the Eastern one (Deuss et al., 
2010). 

In the 80's and 90's, anisotropy was thought to be the 
prominent feature. This is commonly attributed to the 
preferred orientation of iron crystals, possibly acquired 
during solidification but most probably resulting from creep 
flow. Thermal convection, developing a flow characterized 
by a spherical harmonic degree I — 1 (a Clementine shape) , 
orientated along the spin axis of the Earth, was one of 



the first candidates to account for anisotropy (Jeanloz and 
Wcnk, 1988). In this model and the following ones (Weber 
and Machetel, 1992; Buffett, 2009), the nature of the inner 
core boundary (i.e. a phase change) was not considered, a 
classical impermeable boundary condition being preferred. 
With permeable boundary conditions — a phase change 
does not prevent material transfer — , the expression of 
convection at the first harmonic degree is a constant ve- 
locity field across the inner core, that obviously does not 
produce any deformation, but implies melting on one side 
and crystallization on the opposite side. This peculiar sit- 
uation, called translation of the inner core, has recently 
been put forward to explain the hemispherical asymmetry 
of velocity and attenuation (Monnereau et al., 2010). 

If the inner core grows in a superadiabatic regime, 
which is the condition for the onset of convection, an un- 
stable thermal stratification develops, so that any infinites- 
imal thermal heterogeneity of harmonic degree / = 1 (i.e. 
one side colder than the other one) will be amplified. Such 
a heterogeneity will induce a displacement of the inner 
core to maintain its center of mass — shifted toward the 
denser/colder hemisphere — at the center of the Earth. 
The inner core acquires a positive topography on the hot- 
ter and lighter side and a negative one on the opposite side. 
This topography is thermodynamically unstable: the side 
emerging above the phase change melts, bringing hotter 
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material up to the surface, while the sinking side allows 
crystallization. The phase change acts to remove the to- 
pography, which is continuously rebuilt by isostatic equi- 
librium. This feedback results in a permanent drift from 
the crystallizing side to the melting side, the drift velocity 
being controlled by the ability of the outer core to restore 
the adiabatic condition at the surface of the inner core. 

Translation is consistent with multiple scattering mod- 
els of wave propagation. If they do not experience defor- 
mation, iron crystals grow as they transit from one hemi- 
sphere to the other. Larger crystals constituting a faster 
and more attenuating medium, a translation velocity of 
some cm/yr (about ten times the growth rate) is enough 
to account for the superficial asymmetry observed for P- 
wave velocity and attenuation, with grains of a few hun- 
dred meters on the crystallizing side (West) growing up to 
a few kilometers before melting on the East side, and a 
drift direction located in the equatorial plane. 

Translation was also proposed to be responsible for the 
formation of a dense layer at the bottom of the outer core, 
since the high rate of melting and crystallization would 
release a liquid depleted in light elements at the surface 
of the inner core (Alboussiere et al., 2010). This would 
explain the anomalously low gradient of P wave velocity 
in the lowermost 200 km of the outer core (Poupinet et al. , 
1983). 

Clearly, translation cannot account for anisotropy. How- 
ever, convective modes developing at higher harmonic de- 
grees (I > 1) will necessarily induce shear deformation that 
could be a source of anisotropy. The development of these 
modes depends on the Rayleigh number, that controls the 
vigor of the convection, and thus mainly on the viscosity 
of the solid portion of the core. At high Rayleigh number 
(low viscosity), these modes can be dominant and dissi- 
pate the degree I = 1 of the thermal heterogeneities: the 
source of the translation. Thus a viscosity threshold may 
be expected below which translation would not take place. 
This may constrain the viscosity of iron at the conditions 
of the inner core, based on seismological observation. In 
this paper, we present dynamics model of inner core tak- 
ing into account the phase change boundary and study the 
interaction of translation and convection. 

Model setup 

Inner core dynamics obey similar governing equations 
to those used for mantle convection. The specificity lies 
in the boundary conditions required to treat the phase 
change at its surface. 

Momentum equations 

Inertial forces can be neglected because of the high vis- 
cosity of solid iron at the inner core temperature and pres- 
sure conditions, which is at least 10 16 Pas (Yoshida et al., 
1996). Conservation of momentum just expresses the bal- 
ance between buoyancy forces and viscous dissipation. It 



is time independent. 
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t is the deviatoric stress tensor, p the pressure, p the den- 
sity and g the gravity. The permeable surface condition 
is introduced. It describes the balance between the radial 
stress and buoyancy forces induced by topography: 



P\Ria = (Pi - Ps)gicbh, 



(2) 



where r rr is the deviatoric radial stress, Ri C the inner core 
radius, p s the density of solid iron, pi the density of liquid 
iron, gi C b the gravity at the surface of the inner core and h 
the topography. In addition, the surface is considered as 
tangential stress free: 



TrS,<j>\R i€ — 



(3) 



At the first spherical harmonic degree 1 = 1, the momen- 
tum equation written for the full sphere with the above 
conditions comes down to a simple isostasy equilibrium 
where thermal heterogeneities are balanced by a topogra- 
phy (V-r = in Eq. 1) (Monncreau et al., 2010). This sin- 
gularity requires a particular treatment. The three compo- 
nents of the first harmonic degree of the topography h\ m 
are directly related to the position of the center of mass 
anomaly caused by the presence of the inner core in the 
outer core: 
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with a the thermal expansion coefficient and 6i m (0 the 
radial profile of the temperature heterogeneity at degree 
I = 1 and order m (see Appendix A for more details). 

Surface heat exchange with the outer core 

The topography is thermodynamically unstable and 
eroded at a rate depending on the vigor of convection 
within the outer core. In the mantle, the position of 
mineral phase transitions is mainly related to the ambi- 
ent temperature, latent heat exchanges having almost no 
effect. For instance, a transition with a positive Clapey- 
ron slope like olivine to spinel occurs deeper in ascending 
(hot) currents than in dipping slabs. The reverse situa- 
tion happens for the ICB. The turbulent flow in the liq- 
uid maintains the temperature above the surface close to 
the adiabat so that the topography only depends on the 
temperature variations induced by the latent heat effects 
and not on temperature anomalies within the inner core. 
Topography is thus positive where material exits the in- 
ner core because of the cooling induced by the latent heat 
consumption. The local thermodynamical equilibrium is 
achieved when the rate at which the inner core consumes 
or releases latent heat equals the rate at which the outer 
core brings or takes the energy to maintain the adiabatic 



temperature. This can be written as (Alboussiere et al. 
2010): 



V r (Ric) = Fh, 



where 



-pWicb 



fdT n 



aTo 



F 



V dp Pi C p 



uiC p 



L 



(5) 



(0) 



Vr(Ric) is the radial velocity across the phase change, T m 
the melting temperature of iron, T$ the adiabatic temper- 
ature, ui the amplitude of the outer core convective flow at 
the surface of the inner core, C p the specific heat and L the 
latent heat of the phase change. We studied the kinetics of 
this equilibrium in a simple ID model (see Appendix B) in 
order to check that it can be considered as instantaneous. 
The temperature adjustment is reached in a few thousands 
years which is at least one order of magnitude smaller than 
the characteristic time of the convective fluctuations at the 
highest Rayleigh number we considered. 

Since Eq. 5 is linear and time independent, it applies in 
the spectral domain, so that the convection and isostatic 
problems can be separated. The velocity of translation 
is directly obtained from h\ m (Eq. 4) and topography at 
higher degree (I > 1) is computed through the resolution 
of Eqs. 1-3. 

Energy equation 

The inner core dynamics is a moving boundary prob- 
lem that can be treated in all its complexity (Dcgucn and 
Cardin, 2011). For the sake of simplicity, we may neglect 
the variation of the radius with time, and focus on the 
present time dynamics. This assumption is plainly justi- 
fied since the translation velocity required to account for 
seismic properties of the inner core, but also for the for- 
mation of a dense layer above ICB, should exceed the in- 
ner core growth rate (by one or two orders of magnitude) 
(Monnereau et al., 2010; Alboussiere et al., 2010). The en- 
ergy equation is thus written in terms of the temperature 
relative to the adiabat anchored at the ICB (9 = T — T$), 
in which the decrease of ICB temperature with time plays 
the role of an internal heating: 
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where v is the velocity vector, k the conductivity, $ the 
internal heating rate and A the Laplacian operator. The 
surface is assumed isothermal, 9 = 0. This comes to ne- 
glect the small temperature perturbation due to the latent 
heat (~ 10 _3 K) compared to the one involved in the dy- 
namics (~1K). Indeed, the latent heat effects play a much 
more important role in the momentum equation than in 
the energy equation. 

Parameters 

$, the internal heating, results from a competition be- 
tween the heat lost by conduction along the adiabat and 



the secular cooling of the core: 



$ = kAT s - p 8 C p 



dTs 
dt 



(8) 



The superadiabatic regime, required for the onset of 
translation or convection, is defined by $ > 0. Uncertain- 
ties on both the conductivity and the cooling rate con- 
tribute to the uncertainty on $. The former was evalu- 
ated around 60 W/mK (Stacey and Anderson, 2001) ; an 
estimate that has been recently revisited and reduced by a 
factor of 2: 36 W/mK (Stacey and Loper, 2007). dT s /dt 
is proportional to the inner core growth rate Ri c : 
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Ri C being itself proportional to the heat flux at the core- 
mantle boundary (CMB) (Labrossc, 2003), whith walues 
for the latter estimated to be 6 TW and 14 TW (Turcottc 
and Schubert, 2002). It results that $ can be as large as 
2.5 x 10" 8 W/m 3 . 

The vigor of convection is controlled by the Rayleigh 
Number: 
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written here for pure internal heating in a spherical system. 
At first sight, the variation of radius appears of primary 
importance for the evolution of the dynamics of the inner 
core along its history, notably because of the successive 
development of convective modes at higher and higher de- 
gree as the inner core grows (Deguen and Cardin, 2011). 
However, as mentioned before, we focus on the present 
day dynamics. In this case, the least constrained parame- 
ters are the internal heating, $ and the viscosity, n, both 
subject to uncertainties of several orders of magnitude. 

A minimum viscosity of 10 16 Pas was deduced from 
attenuation of seismic waves in the inner core and a max- 
imum of 10 21 Pas was inferred from conjectures on iron 
rheology close to its melting point for pressures at the cen- 
ter of the Earth (Yoshida et al., 1996). Since translation 
does not involve deformation, it is not sensitive to 77, but 
only to $. Their effects have to be studied independently. 

The third and last parameter is the amplitude of the 
convective velocity within the outer core, ui, that acts only 
on the translation velocity, the effect of the topography re- 
maining very small on the convective circulation. Other 
parameters in the expression of the factor F (Eq. 5) remain 
much more constrained than ui. The average convection 
velocity in the outer core was assessed to be 10 -4 m/s 
with 10~ 3 m/s maxima, from the secular variations of the 
geomagnetic field (Amit and Olson, 2006). It is possible 
that a dense layer forms at the base of the outer core (Al- 
boussiere et al., 2010; Buffett, 2011), convection velocity 
could be lower at the surface of the inner core. We will 
also test a velocity ten times lower than the average. 



Table 1: Value of parameters used in the numerical experiments. 



Symbol 


Name 


value 


unit 


a 




Thermal expansion coefficient 


2.5 x 1(T 5 


K" 1 


L ■■ i> 




Heat capacity 


800 


J/Kkg 


T s 




Adiabatic temperature anchored at ICB 


5500 


K 


k 




thermal conductivity 


34 


W/K m 


8T 


J dp 


melting curve slope 




K/Pa 


L 




Latent heat 


800 


J/kg 


Ps 




solid iron density 


11800 


kg/m 3 


Ps 




liquid iron density 


11200 


kg/m 3 


**>ic 




inner core radius 


1220 


km 


9icb 




gravity at the ICB 


4 


m/s 2 


u, 




Outer core velocity 


10~ 5 - 10~ 3 


m/s 


$ 




Internal heating 


2.5 x 10~ 10 - 2.5 x 10~ 8 


W/m 3 


'/ 




Viscosity 


10 16 - 10 20 


Pas 


Ra 




Rayleigh number 


6 x 10 4 - 6 x 10 9 





Numerical methods 

The momentum equation is solved using a classical 
spherical harmonics expansion together with a radial hnitc 
difference solver for each harmonic. The conservation of 
energy is solved with a second order finite volume method. 
More details on the code may be found in Monnereau and 
Yuen (2002). For all calculations we have employed 128 
spherical harmonics and 128 points in the radial direction. 

Results 

Table 1 summarizes the range of investigation for the 
parameters $, 77 and U[. Figure 1 presents the velocity 
fields of experiments with an intermediate value for the 
outer core velocity, m = 10~ 4 m/s. A rapid overview 
shows that translation dominates at high viscosity and 
high internal heating (Fig. lc, e-i), and is absent from ex- 
periments at low viscosity and low internal heating (Fig. la, 
b & d) . Translation is characterized by a surface radial 
velocity, i.e. the radial velocity through the ICB, positive 
on one hemisphere and negative on the opposite, inducing 
melting and crystallization respectively. 

Viscosity effect 

At high viscosity, r\ — 10 20 Pas, translation seems to 
be the sole active mode of heat transfer. Nevertheless, the 
velocity field does not coincide with a pure solid motion, 
except at very low internal heating rate (Fig. lg). It also 
contains a weak degree I = 2 component, attesting the de- 
velopment of convective modes at degree I > 1. This fea- 
ture is not visible on the surface radial velocity maps, but 
is revealed by the non-collinearity of the velocity vectors 
depicted on the cross section of Fig. lh. At lower viscosity, 
r) = 10 18 , translation does not develop for the lowest inter- 
nal heating rate Fig. Id. A degree I = 2 dominates in this 
case, where a sheet-like downwclling draws a great circle. 



When translation still takes place (Fig. le & f), it is asso- 
ciated with one or several planar dipping currents. These 
short wavelengths structures only form within the crystal- 
lizing hemisphere. At low viscosity, 77 = 10 16 Pas, these 2D 
currents give way to cold spots. Translation disappears at 
low and intermediate internal heating rates (Fig. la & b). 
It persists at the highest heating rate, but, contrary to the 
situations at 77 = 10 18 Pas, high degree convective struc- 
tures are present on both crystallizing and melting hemi- 
spheres even if an asymmetry in the number of structures 
is still present. 

In summary, a low viscosity enhances the vigor of the 
convection and favors the development of short wavelength 
structures, a result expected since it corresponds to an 
increase in Rayleigh number. 

Internal heating effect 

Internal heating also enhances convection, as evidenced 
by the increase of the number of structures. However, 
when degree I = 1 dominates the circulation, the major 
effect is an amplification of the translation, whose velocity 
seems only dependent on the internal heating rate. For in- 
stance, the comparison of situations depicted in Fig. le & h 
and Fig. If & i. shows that translation reaches the same 
velocity at r] = 10 20 Pas and 77 = 10 18 Pas, irrespective of 
the convective pattern. This highlights the fact that the 
Rayleigh number is not an independent parameter of the 
coupled system. Besides, this system expresses differently 
at equal Rayleigh numbers as in Fig. 2a where point like 
cold currents inhibit the translation and in Fig. 2f where 
the reverse prevails. Internal heating has a promoting ef- 
fect on the translation. 

Translation vs convection 

Translation induces an asymmetry of the thermal bound- 
ary layer with a thickening on the crystallizing side and 
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Figure 1: Effects of viscosity and internal heating on the style of convection in the inner core. Surface radial velocity fields (hammer projection) 
are displayed for experiments with 3 values of viscosity r) and internal heating <t>. They are snapshots taken after that convection equations 
were integrated to a quasi steady-state regimes. Equatorial cross-sections of the velocity field are also shown for the intermediate internal 
heating value <I> = 2.5 X 10 — 9 Wm — 3 . The background color is function of the velocity amplitude. Note the predominant influence of viscosity, 
with complex convection modes at low viscosity, simple translation mode at high viscosity and hybrid circulation pattern for intermediate 
viscosity. The Rayleigh number, Ra, is indicated in blue to the upper left corner of each plot. Identical Rayleigh numbers, corresponding 
here to identical <J>/»? ratio, may result in different convective pattern, because internal heating has a larger impact on translation than on 
convection. 



a thinning on the melting side. When this thinning is 
enough to reduce the local Rayleigh number of the ther- 
mal boundary layer below the critical value, no convective 
instability can begin. This leads to situations where con- 
vective structures are restricted to the crystallizing side as 
in Fig. le & f. This also explains the large difference be- 
tween Fig. Id & i, both corresponding to experiments with 
a Rayleigh number more than 200 times the critical value 
for the onset of convective structure at a spherical har- 
monic degree I = 2 (2607 with impermeable surface con- 
ditions (Chandrasckhar, 1961)). At high Rayleigh num- 
ber, when the conditions of this thinning are not fulfilled, 
higher harmonic degree structures, such as coldspots, de- 
velop. These uniformly distributed downwellings suppress 
the thermal anomaly at harmonic degree I = 1 that trans- 
lation requires to develop. Thus a real competition is play- 
ing between convection and translation. 

These are two distinct modes of heat transfer, which 
can be verified by examining global quantities. For in- 
stance, horizontally averaged temperature profiles exhibit 
very different shapes depending on whether translation op- 



erates or not. The classical temperature profile obtained 
for convection in a full sphere is marked by a top thermal 
boundary layer and a maximum located just below. When 
the translation develops, this shape evolves towards an el- 
liptic integral profile E(tt, r) where the maximum is now 
reached at the center (see in Fig. 2 changes from profiles b 
to e, then f). This is the analytical solution of the temper- 
ature equation when conduction and time dependence are 
neglected. An additional major effect is the strong cooling 
produced by the translation as illustrated by the compar- 
ison of profiles of experiments run at the same Rayleigh 
number (a and f, or d and i). The permeable boundary 
that allows advective heat transfer across the surface is re- 
sponsible for this cooling. This result was already shown in 
a former study, where such permeable conditions were used 
to mimic the cooling due to intense volcanism on Io (Mon- 
nereau and Dubuffct, 2002). In this case, the exponent of 
the power scaling relation between averaged temperature 
and Rayleigh number is —1/2 with permeable conditions, 
instead of the classical value —1/4 found with impermeable 
conditions (McKcnzic ct al., 1974). 



The averaged temperature decreases proportionally with 
Ra~ 1 ' 2 only when translation dominates the circulations 
(see Fig. 3). Elsewhere the power low exponent is —0.2. 
Actually, the non dimensional averaged temperature is 
just proportional to $~ 1//2 , since viscosity has no effect 
on translation. This relation can be easily understood by 
considering that, when translation dominates, the energy 
equation comes down to a balance between internal heat- 
ing and advection, so that the product of velocity and 
temperature is proportional to $. This also expresses the 
fact that all the energy produced internally is evacuated 
by advection. As the translation velocity is proportional 
to the temperature (see Eqs [4] & [5]), the temperature 
is proportional to $ x / 2 . Finally, the non dimensional tem- 
perature is obtained by dividing the temperature by the 
temperature scale <&i? 2 c /6fc. When the geometry does not 
allow the translation mode, as in the case of convection 
between spherical shells, or for all convective modes at de- 
gree I > 1, the velocity is proportional to $/r/, so that 
the non-dimensional temperature varies in proportion to 
Ra- 1 ' 2 . 

The origin of the difference between these modes lies 
in the fact that the permeable boundary conditions act 
differently at degree I — 1 and degrees I > 1. In the case 
of translation, the radial velocity across the ICB is simply 
proportional to the topography through Eq. [5]. For con- 
vection (I > 1), the topography is also proportional to the 
vertical stress Eq. [2] and thus to dv r /dr. This last rela- 
tion expresses the reverse of Eq. [5] : in absence of phase 
change, i.e. when uiO, the topography is maximal when 
the velocity across the ICB is zero. The antagonism of 
both relations reduces the amplitude of the radial velocity 
at degree I > 1. Further, since the stress is proportional 
to the viscosity, the ICB becomes impermeable at low vis- 
cosity, as observed in Fig. 4. 

Temporal evolution 

Most of the circulations shown in Fig. 1 are time de- 
pendent except those run with a viscosity r\ ^ 10 20 Pas. In 
a relatively limited range of viscosities, we observe periodic 
regimes. For instance at 10 19 Pas the energy is exchanged 
between degree 2 to degree 1 in an oscillating regime of pe- 
riod of 1.5 Gyr. Such a period comparable with the age of 
the inner core age is irrelevant because inner core growth 
is not modeled. At 10 18 Pas translation evolves pseudo- 
steadily while higher degrees develop chaotically. Circula- 
tion becomes completely chaotic below 10 18 Pas. The time 
scale of variation are, in these cases, smaller than 100 Myr 
and thus appears rapid at the scale of the inner core age. 

Effect of convection velocity in the outer core 

Convection velocity ui drives material transfer through 
the phase change, owing to its capacity to extract or bring 
latent heat. It consequently drives surface radial velocity 
and mainly affects the translation. For instance at low 
viscosity (r) = 10 16 Pas), we observed that increasing ui 
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Figure 2: Non dimensional temperature profile of experiments shown 
in Fig. 1. Labels are the same as in Fig. 1. The non dimensional 
temperature can be scaled by $i?? c /6fc, where the conductivity, k, 
has been set to 36 W/mK. The line style is related to the internal 
heating rate: dotted lines for $ = 2.5 X 10~ 10 W/m 3 , dot-dashed 
lines for $ = 2.5 X l()- 9 W/m 3 , plain lines for <J> = 2.5 X 10~ 8 W/m 3 . 
The color changes from dark blue to red as the Rayleigh number 
increases. Rayleigh numbers are indicated on Fig. 1. For instance, 
profiles d and i, or a and f correspond to experiments with the same 
Rayleigh number, 6 X 10 s and 6 X 10 7 respectively. A comparison 
of these profiles shows the strong cooling induced by translation. 
Profiles resulting from almost pure translation are parabolic (i, h & 
g) and evolve toward the classical shape observed for convection in 
a full sphere (a and b). 



has almost no impact on convection, except on the surface 
radial velocities, which are proportional to ui. But neither 
the shape or the number of structures are affected. On 
the other hand, its enhancing effect on the translation is 
clear. This one is absent at ui = 10~ 5 m/s, just appears at 
Ui = 10 _4 m/s for an internal heating rate $ = 2.5 x 10 -8 
W/m 3 , with convective structure equally distributed on 
both hemisphere (Fig. lc). At ui = 10~ 3 m/s, translation 
is present whatever the internal heating rate, but domi- 
nates for an internal heating rate $ ^ 2.5 x 10 -9 W/m 3 , 
as revealed by the asymmetric repartition of convective 
structures, 

Fig. 5 summarizes experiments performed in this study. 
Globally, it highlights the fact that i) viscosity weakens the 
convection, ii) the liquid velocity at the inner core surface 
promotes translation and Hi) internal heating promotes 
both convection and translation, but with a larger effect 
on the latter. We can see on this figure that most of the 
situations found at given values of 77, $ and ui are equiv- 
alent to the ones found at 77/ 10, $/10, and lOu; 

Implication for the viscosity of the inner core 

The seismic asymmetry of the inner core and particu- 
larly the correlation of velocity with attenuation of com- 
pressive wave traveling in the uppermost part of the inner 
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Figure 3: Non dimensional averaged temperature as a function of 
the Rayleigh number of experiments shown in Fig. 1. Labels are the 
same as in Fig. 1. 
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Figure 4: Effect of viscosity on surface velocities (left) and maximum 
velocities (right). A threshold around 10 19 Pas is clear, with constant 
velocities above denoting that solid translation prevails. Below, the 
surface radial velocity decreases whereas the surface horizontal ve- 
locity increases, denoting a change of the ICB from permeable to 
non-permeable state. The maximum velocity reveals also this change 
showing a constant exchange with the outer core above 10 19 Pas, and 
a power law dependence characteristic impermeable conditions below 
(the slope is around -2/5). These experiments have been run with 
U { = lCT 4 m/s and $ = 2.5 X 10~ 9 W/m 3 



core has been interpreted as the signature of iron crystal 
growth during their transit from a freezing side (facing 
Peru) to a melting side (facing Indonesia). This model re- 
quires that translation be dynamically possible. A first 
condition is that the inner core be in a superadiabatic 
state, corresponding here to $ > (see Eq. [8]). A second 
condition concerns the viscosity of solid iron at inner core 
conditions, 77. If the viscosity is too low, high harmonic 
degree convection takes place, preventing the development 
of the translation. Whatever the internal heating consid- 
ered here (i.e. the core cooling rate), almost pure trans- 
lation prevails, above a viscosity threshold that only de- 
pends on the vigor of the outer core convection. In Fig. 5, 
this threshold appears where the constant product of the 
average liquid iron velocity above the ICB ui and the vis- 
cosity exceeds 10 16 Pam. Since the most accepted value 
for ui lies around 10 m/s, this threshold corresponds to 
i] = 10 20 Pas. 



16 

17 

. 18 

19 

I 

20 
21 



© © © 16 

© © © 17 

• © • 18 

© © © 19 

O O <D 20 

€ C € 



21 l 



1 

- © 


1 

© 


1 

© 


- © 


© 


© 


- • 


© 


© 


- • 


© 


© 


- • 





© 


L € 





© 





0.25 


2.5 


25 


16 


1 

- © 


1 

© 


1 

© 


17 


e 


© 


© 


18 


- © 


© 


© 


19 


- • 


© 


© 


20 


- © 


© 


© 


21 


L € 


© 


O 



ul = 10 m/s 



ul = 10 m/s 



ul = 10 m/s 



^B 1-1 dominates (steady state) 1-1 (pseudo-steady) + high degrees ^p 

on one hemisphere (chaotic) 
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^p 1=1, 1=2 comparable amplitudes 
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Figure 5: Summary of the convection patterns obtained in this study. 



Just below, at rj = 10 19 Pas, harmonic degree I = 2 
establishes or alternates with the harmonic degree I = 1 
(translation) in a periodic or quasi periodic regime. Be- 
cause of their time scale exceeding the age of the inner 
core, these situations appear irrelevant and would necessi- 
tate a modeling that takes into account inner core growth. 
In any case, a circulation dominated by a degree / = 2 
would not satisfy the observed seismic asymmetry. 

Translation is recovered at viscosity one order of mag- 
nitude lower, 77 = 10 18 Pas, but embedded with short 
wavelength convective structures that develop on the freez- 
ing side only (see cross section of Fig. le). The resulting 
deformation precludes crystal growth, but may imprint a 
texture that could be responsible for the seismic anisotropy 
of the inner core. Seismic anisotropy may be subject to a 
hemispherical variation, with the western hemisphere dis- 
playing much stronger anisotropy than the eastern hemi- 
sphere (Tanaka and Hamaguchi, 1997; Deuss et al., 2010). 
On Fig. le, the hotter hemisphere is not sheared by convec- 
tive currents, material being simply advected toward the 
melting side. Annealing during this transit would erase the 
texture and allow crystal growth (Bergman ct al., 2010). 
This asymmetric convection is still present at lower vis- 
cosity, 77 ^ 10 17 Pas, but only for values of ui ten times 
larger than the classically expected at the base of the outer 
core. Such a vigorous circulation appears in opposition 
with the possible presence of a stable high density layer 
above the ICB, produced by the fusion of inner core ma- 
terial depleted in light elements, for which the translation 
of the inner core has also been invoked (Alboussiere et al., 
2010). For the expected ui value, a degree I = 1 remains 
for 77 ^ 10 17 Pas, but short wavelength convective struc- 
tures stir both hemispheres preventing any asymmetry of 



crystal size and orientation and the formation of a dense 
layer at the surface of the inner core. 

In this study, the external environment was considered 
as homogeneous. Global circulation in the liquid core may 
induce variations of the thermodynamic conditions at the 
surface of the inner core, which are able to significantly 
disturb its growth rate (Aubert et al., 2008; Gubbins et al., 
2011). External forcing (Yoshida et al., 1996; Aubert et al., 
2008; Dcgucn et al., 2011) may act in conjunction with 
thermal convection and may stabilize degree 2 and degree 
1 together. 

If translation generates the seismic asymmetry of the 
inner core and if its anisotropy results from convection, 
our numerical simulations show that viscosity is comprised 
between 10 18 and 10 20 Pas for reasonable values of con- 
ductivity and age of the inner core. For viscosities lower 
than 10 18 Pas the "sluggish" inner core is impermeable, 
meaning that exchanges with the outer core are small, it 
is deformed and has no hemispherical pattern. For viscosi- 
ties higher than 10 20 Pas the rigid inner core is permeable, 
non-deformed but hemispherical. For intermediate viscosi- 
ties (10 18 Pas) particular conditions implying a boundary 
layer close to its Rayleigh number, allows convection to 
develop over a bulk translation. This configuration pos- 
sibly makes thermal convection a process able to explain 
seismic dichotomy and large-scale variation of anisotropy. 

Appendix A. Force balance of the inner core 

Convection within the inner core affects the density 
field. The harmonic degree I = 1 of these variations will 
shift the center of mass of the inner core. To visualize 
this shift, let us consider the inner core as a solid sphere 
of average density p s , entirely covered with a liquid ocean 
(the outer core) of density pi (see Fig. A. 6). The density 
of the solid, marked by a variation at a spherical harmonic 
degree I = 1 due to a temperature variation, is: 



i-a J2 @ r(r)ir(0,< 



(A.l) 



where @™(r) are the components of degree I = 1 of the 
temperature field. More simply, we may write: 



p = p s [l -~aOi(r)cos9], 



(A.2) 



where 9 is the angle made by the radius originating at 
the center of the figure of the solid C and the direction 
of the axis defined by C and the center of mass of the 
system O, 6i(r) being the radial profile of the temperature 
heterogeneity at degree 1 = 1. The ocean surface is an 
equipotential of the gravity field and thus is centered on 
the center of mass O. The position of the center of mass 
is the one that reduces to zero the first torque of the mass 
anomaly that represents the inner core within the outer 




Figure A. 6: 



core, the contribution of the outer core being zero: 

/ OMp(M)dv = / OMpidv + 
Jv Jv 

/ OM [ Ps - ap a @i(M) - pi] dv = 

./inner core 



(A.3) 



Introducing the radial vector OM = OC + r leads to: 



2ir ix Rio 



(OC + r cos 9) [p s — pi — p s aQi(r) cos 9] 







(A.4) 



and finally: 



OC 



p s a 



RUPs -Pi) Jo 



r 2 sm9drd9d(/) = 0, 



r 3 6i(r)dr. 



(A.5) 



The distance OC also corresponds to the amplitude hi of 
the harmonic degree I = 1 of the topography of the solid 
sphere referenced at O ; this allows to relate /i™ to 0" 1 
through the coordinate of O relatively to C. 

Appendix B. Phase change kinetics ID model 

A delicate point in solving this problem concerns the 
computation of the inner core topography. At the har- 
monic degree 1 = 1, the topography is given by the offset 
position of center of mass, whereas it results from the con- 
vection velocity, at degrees I > 1. The solution would 
consist in writing an integral boundary condition via a to- 
pography differential equation: 



dh 



Vtf (R ic ) ■ Vh = V r {Ric) + Vf, 



where vh(Ric is the surface horizontal velocity, v r (Ri C ) 
the radial velocity across the ICB and vj the solidifica- 



tion/melting velocity. If the kinetics of the phase change 
is fast, the above equation comes down to a simple bal- 
ance between v r (Ri C ) and Vf, so that the topography can 
be calculated in the spectral domain. In order to ver- 
ify that adjustment of melting front is almost instanta- 
neous compared to inner core convection time scale, we 
developed a ID model of the time evolution of the melt- 
ing/solidification front. 

Appendix B.l. Model 

We solve the energy equation in which the latent heat, 
L, has been introduced as in Christensen and Yuen (1985) 




DT DT 

P C P —- P L—+aTpgv z 



, d 2 T a . , 



(B.l) 



r is a function describing the solid fraction: 



1 



r ( 7r ) = o(l + tanh(7r)), with n 



p-Po-jT 
Sp 



(B.2) 

n is the pressure offset due the phase change, normalized 
by its pressure thickness dp. 7 is the Clapeyron slope of 
the pause change, p the pressure, po a reference pressure. 
S(z), the source/sink term, corresponds to the capacity of 
outer core to provide or extract the phase change latent 
heat: 

JT"! 

S(z) = -pC p (T - T ad )ui- 



with 



S becomes: 
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Expanding DT/Dt: 
DT 

rn 



dT ( dn DT 
dn' {dfrJt 



dn Dp 
dp~Dt 



(B.3) 
(B.4) 

(B.5) 
(B.6) 



and introducing the temperature dependence of latent heat: 



L = 



P z 



(B.7) 



allows to write the specific heat C' and the thermal ex- 
pansion coefficient a' modified by latent heat (Christensen 
and Yuen, 1985): 

r»(\ n . T(z) 1 2 5pdT . IP dT 

C p (z) = C p + — — et a'{z)=u+— — . 

y p z op an pop an 



Finally conservation of heat (B.l) becomes: 

pc' p 



dT dT 

dt dz 



dz 2 



+ a'pgv z T = k 1 ^ + S{z). (B.8) 



Figure B.7: Temperature profile evolution for simulation 1. Initial 
temperature profile is in blue thick line. Iron melting curve is in 
green thick line. Temperature profile at thermal equilibrium is red. 
Shaded orange tons correspond to intermediate profiles. 



We use a Backward Euler scheme to solve differential 
equation (B.8). Boundary conditions arc differently de- 
fined when material enters or leave the box. Solid Iron 
(the inner core) is supposed to occupy the bottom half 
of the domain. In the case of solidification (v z < 0), the 
top temperature is ascribed to the outer core temperature, 
which is T top = T a d- At the bottom, we impose the tem- 
perature gradient to be equal to the adiabatic gradient. In 
the case of melting (v z > 0), the bottom temperature is 
set to the adiabatic temperature plus an anomaly resulting 
from internal heating: T top — T aa - + ST. At the top, we 
impose the temperature gradient to be equal the adiabatic 
gradient. 

Results 

We performed five simulations varying the vertical ve- 
locities v z -- supposed to be the radial velocity across 
the Inner core Bounadry (ICB) — , the convective velocity 
within the liquid outer core just above the ICB and the 
temperature within the inner core. Table B.2 summarizes 
the characteristics of these experiments. 

Fig. B.7 show the time evolution of temperature pro- 
files in case of melting by decompression (v z > 0). We 
see that the consumption of latent heat cools the inner 
core by a few hundredths of kelvin, shifting up the phase 



Table B.2: Numerical experiment characteristics. v z is the vertical 
velocity, n; convection velocity inthe outer core, h the topography at 
thermal equilibrium. 

Exp. v z (cm/yr) ui (m/s) Initial T(z) (h) (m) 
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Figure B.8: Topography evolution for the five simulations. 

change by a hundred meters. This state is reached in less 
that 10 4 years, which is shorter than time scale involved in 
convective process. A symmetric situation is observed in 
the case of solidification when v z < 0. The time necessary 
to reach a steady state depends on the difference between 
the temperature of the inner core material and the adia- 
bat. This time is four time larger for IK of difference. On 
the other hand, we find, as expected, that the topography 
reached is non dependent on this parameter and that it is 
just proportional to v z and to 1/ui 

Equilibrium specific time, on the order of 10 4 to 10 5 
years, remains short enough when compared to convection 
time scales to consider thermal equilibrium of the topog- 
raphy as instantaneous. 
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